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Abstract 

Charm production at the LHC in pp collisions at i/s = 7TeV is studied with the 
LHCb detector. The decays Z?°— » K~w + , D + -> K-^+n+, L>*+ -> D°(K-7r+)7r+, 
— > 4>{K~ K + )tt + , 4+— » p-ftT _ 7r + , and their charge conjugates are analysed in a 
data set corresponding to an integrated luminosity of 15 nb -1 . Differential cross- 
sections da/dpT are measured for prompt production of the five charmed hadron 
species in bins of transverse momentum and rapidity in the region < p^ < 8 GeV/ c 
and 2.0 < y < 4.5. Theoretical predictions are compared to the measured differential 
cross-sections. The integrated cross-sections of the charm hadrons are computed 
in the above pr-y range, and their ratios are reported. A combination of the five 
integrated cross-section measurements gives 

o"(cc) PT <8G e v/c,2.o< ?; <4.5 = 1419 ± 12 (stat) ± 116 (syst) ± 65 (frag) \ib, 

where the uncertainties are statistical, systematic, and due to the fragmentation 
functions. 
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1 Introduction 



Measurements of the production cross-sections of charmed hadrons test the predictions of 
quantum chromodynamic (QCD) fragmentation and hadronisation models. Perturbative 
calculations of charmed hadron production cross-sections at next-to-leading order using 
the Generalized Mass Variable Flavour Number Scheme (GMVFNS) [j}j6] and at fixed 
order with next-to-leading-log resummation (FONLL) [7}{To] reproduce the cross-sections 
measured in the central rapidity region (|y| < 1) in pp collisions at \/s = 1.97 TeV at the 
Fermilab Tevatron collider 11 and the cross-sections measured in the central rapidity 



region (\y\ < 0.5) in pp collisions at y/s = 2.96 TeV [12] and at yfs = 7 TeV 13 14 
at the CERN Large Hadron Collider (LHC). The LHCb detector at the LHC provides 
unique access to the forward rapidity region at these energies with a detector that is 
tailored for flavour physics. This paper presents measurements with the LHCb detector 
of D°, D + , D+, D* + , and A£ production in the forward rapidity region 2.0 < y < 4.5 in 
pp collisions at a centre-of-mass energy of 7 TeV. Throughout this article, references to 
specific decay modes or specific charmed hadrons also imply the charge conjugate mode. 
The measurements are based on 15 nb -1 of pp collisions recorded with the LHCb detector 
in 2010 with approximately 1.1 visible interactions per triggered bunch crossing. 

Charmed hadrons may be produced at the pp collision point either directly or as 
feed-down from the instantaneous decays of excited charm resonances. They may also be 
produced in decays of 6-hadrons. In this paper, the first two sources (direct production 
and feed-down) are referred to as prompt. Charmed particles from 6-hadron decays are 
called secondary charmed hadrons. The measurements described here are the production 
cross-sections of prompt charmed hadrons. Secondary charmed hadrons are treated as 
backgrounds. No attempt is made to distinguish between the two sources of prompt 
charmed hadrons. 



2 Experimental conditions 

The LHCb detector |15| is a single-arm forward spectrometer covering the pseudorapidity 
range 2 < rj < 5, designed for the study of particles containing b or c quarks. The detector 
includes a high precision tracking system consisting of a silicon-strip vertex detector 
surrounding the pp interaction region, a large-area silicon-strip detector located upstream 
of a dipole magnet with a bending power of about 4 Tm, and three stations of silicon-strip 
detectors and straw drift-tubes placed downstream. The combined tracking system has a 
momentum resolution (Ap/p) that varies from 0.4% at 5GeV/c to 0.6% at lOOGeV/c and 
an impact parameter (IP) resolution of 20 \xm for tracks with high transverse momentum. 
Charged hadrons are identified using two ring-imaging Cherenkov detectors. Photon, 
electron, and hadron candidates are identified by a calorimeter system consisting of 
scintillating-pad and pre-shower detectors, an electromagnetic calorimeter, and a hadronic 
calorimeter. Muons are identified by a system composed of alternating layers of iron and 
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multiwire proportional chambers. The trigger consists of a hardware stage, based on 
information from the calorimeter and muon systems, followed by a software stage that 
applies a full event reconstruction. 

During the considered data taking period, the rate of bunch crossings at the LHCb 
interaction point was sufficiently small that the software stage of the trigger could process 
all bunch crossings. Candidate events passed through the hardware stage of the trigger 
without filtering. The software stage of the trigger accepted bunch crossings for which 
at least one track was reconstructed in either the silicon-strip vertex detector or the 
downstream tracking stations. The sample is divided into two periods of data collection. 
In the first 1.9 ±0.1 nb _1 all bunch crossings satisfying these criteria were retained. In the 
subsequent 13.1 ± 0.5 nb -1 the trigger retention rate was limited to a randomly selected 
(24.0 ± 0.2)% of all bunch crossings. 

For simulated events, pp collisions are generated using Pythia 6.4 16 



LHCb configuration 17 that employs the CTEQ6L1 parton densities 



with a specific 
18 . Decays of 

hadronic particles are described by EvtGen [19] in which final state radiation is generated 
using Photos |20|. The interaction of the generated particles with the detector and its 



response are implemented using the Geant4 toolkit 21 as described in Ref. 22 



3 Analysis strategy 

The analysis is based on fully reconstructed decays of charmed hadrons in the following 
decay modes: D ^ R-tt + , D + -> K-ir + n + , D* + -> D°(K-it + )it + , D+ -> <p(K- K + )ti + , 
and A£—} pK~7i + . Formally, the D°— >■ K~n + sample contains the sum of the Cabibbo- 
favoured decays D° — > K~n + and the doubly Cabibbo-suppressed decays D°—> K~n + . 
For simplicity, we will refer to the combined sample by its dominant component. 

The measurements are performed in two-dimensional bins of the transverse momentum 
(pr) and rapidity (y) of the reconstructed hadrons, measured with respect to the beam 
axis in the pp centre-of-mass (CM) frame. For the D°, D + , D* + , and measurements, 
we use eight bins of uniform width in the range < pt < 8 GeV/ c and five bins of uniform 
width in the range 2.0 < y < 4.5. For the A£ measurement, we partition the data in two 
ways: six uniform p^ bins in 2 < p^ < 8 GeV/c with a single 2.0 < y < 4.5 bin and a single 
2 < < 8 GeV/c bin with five uniform y bins in 2.0 < y < 4.5. 

3.1 Selection criteria 

The selection criteria were tuned independently for each decay. The same selection criteria 
are used for D° — > K~n + candidates in the D° and D* + cross-section measurements. We 
use only events that have at least one reconstructed primary interaction vertex (PV). Each 
final state kaon, pion, or proton candidate used in the reconstruction of a D°, D + , D+, or 
A£ candidate must be positively identified. Because of the relatively long lifetimes of the 
D°, D + , D+, and hadrons, the trajectories of their decay products will not, in general, 
point directly back to the PV at which the charmed hadron was produced. To exploit 
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this feature, the selections for these decays require that each final state candidate has a 
minimum impact parameter \ 2 (IP X 2 ) with respect to the PV. The IP \ 2 is defined as 
the difference between the \ 2 °f the PV reconstructed with and without the considered 
particle. For the D° and A£ reconstruction, a common IP \ 2 requirement is imposed on all 
final state particles. For the D + and candidates, progressively stricter limits are used 
for the three daughters. Final-state decay products of charmed hadrons have transverse 
momenta that are generally larger than those of stable charged particles produced at the 
PV. Applying lower limits on the of the final state tracks suppresses combinatorial 
backgrounds in the selections of D°, D + , and A£ samples. 

The selections of candidate charmed hadron decays are further refined by study- 
ing properties of the combinations of the selected final state particles. Candidate 

— > <ft(K~ K + )ir + decays are required to have a K~K + invariant mass within ±20 MeV/c 2 



of the 0(1020) mass 23 . The decay products for each candidate charmed hadron must be 
consistent with originating from a common vertex with a good quality fit. The significant 
lifetimes of D°, D + , D+, and A£ hadrons are exploited by requiring that the fitted decay 
vertexes are significantly displaced from the PV. The trajectory of a prompt charmed 
hadron should point back to the PV in which it was produced. For D° candidates this 
is exploited as a requirement that IP \ 2 < 100. For D° decays, we use one additional 
discriminating variable: the angle between the momentum of the D° candidate in the 
laboratory frame and the momentum of the pion candidate from its decay evaluated in 
the D° rest frame. The cosine of this angle has a flat distribution for D° decays but peaks 
strongly in the forward direction for combinatorial backgrounds. Candidate D* + decays 
are reconstructed from D° and slow pion candidates. Figures [Ifj3] show the invariant mass 
distributions and the log 10 (IPx 2 ) distributions of the selected charmed hadron candidates. 

We factorise the efficiencies for reconstructing and selecting signal decays into com- 
ponents that are measured with independent studies. The particle identification (PID) 
efficiencies for pions, kaons, and protons are measured in data in bins of track p^ and 
pseudorapidity, r], using high purity samples of pions, kaons, and protons from K®, 0(1020), 
and A decays. The effective total PID efficiency for each (ppj y) bin of each charmed 
hadron decay mode is determined by calculating the average efficiency over the bin using 
these final state PID efficiencies and the final state (pt,v) distributions from simulated 
decays. The efficiencies of the remaining selection criteria are determined from studies 
with the full event simulation. 



3.2 Determination of signal yields 

We use multidimensional extended maximum likelihood fits to the mass and log 10 (IPx 2 ) 
distributions to determine the prompt signal yields. For the D* + — > D°tt + mode the 
log 10 (IPx 2 ) of the daughter D° is used. The selected candidates contain secondary back- 
grounds from signal decays produced in decays of 6-hadrons and combinatorial backgrounds. 
The D* + — > D°7i + decay has two additional sources of background from D° decays com- 
bined with unrelated slow pion candidates: prompt random slow pion backgrounds in which 
the D° mesons are produced at the PV and secondary random slow pion backgrounds in 
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Figure 1: Mass and log 10 (IPx 2 ) distributions for selected D° — )> K 7r H 



candidates showing (a) the masses of the D° candidates, (b) the log 10 (IPx 



and D + — > K~it + it + 
2 ) distribution of D 

candidates for a mass window of ±16MeV/c 2 (approximately ±2cr) around the fitted m(K~-K + ) 
peak, (c) the masses of the D + candidates, and (d) the log 10 (IP X 2 ) distribution of D + candidates 
for a mass window of ±HMeV/c 2 (approximately ±2<r) around the fitted m(K~ir + ir + ) peak. 
Projections of likelihood fits to the full data samples are shown with components as indicated in 
the legends. 
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Figure 2: Mass and log 10 (IPx 2 ) distributions for selected D°(K~ir + )ir + candidates 

showing (a) the masses of the D° candidates for a window of ±1.6MeV/c 2 (approximately ±2<r) 
around the fitted Am peak, (b) the differences between the D* + and D° candidate masses for a 
mass window of ±16MeV/c 2 (approximately ±2<r) around the fitted m(K~ir + ) peak, and (c) the 
log 10 (IPx 2 ) distribution of the D° candidate for a mass signal box of ±16MeV/c 2 around the 
fitted m(K~ir + ) peak and ±1.6MeV/c 2 around the fitted Am peak. Projections of a likelihood 
fit to the full data sample are shown with components as indicated in the legend. The l D 
backgrounds' component is the sum of the secondary, prompt random slow pion, and secondary 
random slow pion backgrounds. 
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Figure 3: Mass and l og 10 (IP x 2 ) distributions for selected Df — > <fi( K K + )ir + and A£ ->■ pii" 7r + 
candidates showing (a) the masses of the Df candidates, (b) the log 10 (IPx 2 ) distribution 



of Df candidates for a mass window of ±8MeV/c 2 (approximately ±2<r) around the fitted 
m{(j){K~ K + )tt + ) peak, (c) the masses of the A+ candidates, and (d) the log 10 (IPx 2 ) distribution 
of Af candidates for a mass window of ±12MeV/c 2 (approximately ±2cr) around the fitted 
m(pK~Tr + ) peak. Projections of likelihood fits to the full data samples are shown with components 



as indicated in the legends. 
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which the D° mesons are produced in decays of 6-hadrons. The combinatorial backgrounds 
are separated from the remaining components with the reconstructed D°, D + , Df, and 
mass distributions. Analysis of the log 10 (IP \ 2 ) distributions allow separation of the prompt 
signal and secondary backgrounds. The additional random slow pion backgrounds in the 
D* + — > D (K~7t + )tt + mode are identified in the distribution of the difference Am between 
the masses of the D* + and D° candidates. Thus the prompt signal yields for D°, D + , Df, 
and decays are measured with two-dimensional fits to the mass and log 10 (IPx 2 ), and 
the prompt signal yields for D* + decays are determined with three-dimensional fits to the 
D° candidate mass, Am, and log 10 (IP x 2 )- 

The extended likelihood functions are constructed from multidimensional probability 
density functions (PDFs). For each class of events, the multidimensional PDF is the 
product of an appropriate one-dimensional PDF in each variable: 

Prompt signal: The mass distributions are represented by Crystal Ball functions (24) 
for D° decays (both direct and from D* + mesons), double Gaussian functions for 
the D + and Df modes, and a single Gaussian function for the A£ mode. The 
Am distribution for the D* + mode is represented by a Crystal Ball function. The 
log 10 (IP x 2 ) distributions are represented by bifurcated Gaussian functions with 
exponential tails defined as 

if x < p — pl • a • (1 — e), 
if p — pi ■ cr • (1 — e) < x < p, 
if p < x < p + p R ■ a ■ (1 + e), 
if p + p R ■ a ■ (1 + e) < x, 

(1) 

where p is the mode of the distribution, o is the average of the left and right Gaussian 
widths, e is the asymmetry of the left and right Gaussian widths, and Pl(r) is the 
exponential coefficient for the left (right) tail. 

Secondary backgrounds: The functions representing the mass (and Am) distributions 
are identical to those used for the prompt signal in each case. The log 10 (IPx 2 ) 
distributions are represented by /bg functions. 

Combinatorial backgrounds: The mass distributions are represented by first order poly- 
nomials. The log 10 (IPx 2 ) distributions are represented by /eg functions. The Am 
distribution for the D* + mode is represented by a power-law function C (Am — M 7r ) p 
where the exponent p is a free parameter; is the pion mass and C is a normalisa- 
tion constant. 

Prompt random slow pion backgrounds (D* + only): The functions representing the 
mass and log 10 (IP x 2 ) distributions are identical to those used for the prompt signal. 
The function representing the Am distribution is the same power law function as 
that used for the combinatorial backgrounds. 
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Secondary random slow pion backgrounds (D* + only): The functions representing 
the mass and log 10 (IP x 2 ) distributions are identical to those used for the secondary 
backgrounds. The function representing the Am distribution is the same power law 
function as that used for the combinatorial backgrounds. 

Shape parameters for the log 10 (IPx 2 ) distributions of combinatorial backgrounds are fixed 
based on fits to the mass sidebands. Those of the prompt signal, secondary backgrounds, 
and random slow pion backgrounds are fixed based on fits to simulated events. Figures [l}|3] 
show the results of single fits to the full < pt < 8 GeV/c, 2.0 < y < 4.5 kinematic region. 

The extended maximum likelihood fits are performed for each p^-y bin. We simul- 
taneously fit groups of adjacent bins constraining to the same value several parameters 
that are expected to vary slowly across the kinematic region. The secondary background 
component in the A£ mode is too small to be measured reliably. We set its yield to zero 
when performing the fits and adopt a systematic uncertainty of 3% to account for the 
small potential contamination from secondary production. 

3.3 Systematic uncertainties 

There are three classes of systematic uncertainties: globally correlated sources, sources 
that are correlated between bins but uncorrelated between decay modes, and sources that 
are uncorrelated between bins and decay modes. The globally correlated contributions 
are the uncertainty on the measured luminosity and the uncertainty on the tracking 
efficiency. The former is a uniform 3.5% for each mode. The latter is 3% per final state 
track in the D°, D + , D+, and A£ measurements and 4% for the slow pion in the D* + 
measurement. We adopt the uncertainty of the branching fractions as a bin-correlated 
systematic uncertainty. Systematic uncertainties of the reconstruction and selection 
efficiencies include contributions from the limited size of the simulated samples, failures in 
the association between generated and reconstructed particles in the simulation, differences 
between the observed and simulated distributions of selection variables, and differences 
between the simulated and actual resonance models in the D + and A£ measurements. The 
yield determination includes uncertainties from the fit models, from peaking backgrounds 
due to mis-reconstructed charm cross-feed, and from potential variations in the yields 
of secondary backgrounds. Where possible, the sizes of the systematic uncertainties 
are evaluated independently for each bin. The sources of systematic uncertainties are 
uncorrelated, and the total systematic uncertainty in each bin of each mode is determined 
by adding the systematic uncertainties in quadrature. Table [T] summarises the systematic 
uncertainties. 

As cross-checks, additional cross-section measurements are performed with the decay 
modes D°— > K~7t + 7t~tt + and D + — > <f)(K~ K + )n + and with a selection of D° — > K~n + 
decays that does not use particle identification information. Their results are in agreement 
with the results from our nominal measurements. 
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Table 1: Overview of systematic uncertainties and their values, expressed as relative fractions 
of the cross-section measurements in percent (%). Uncertainties that are computed bin-by-bin 
are expressed as ranges giving the minimum to maximum values of the bin uncertainties. The 
correlated and uncorrelated uncertainties are shown as discussed in the text. 



Source 


D° 


D *+ 


D+ 


Dt 




Selection and reconstruction (correlated) 


1.6 


2.6 


4.3 


5.3 


0.4 


(uncorrelated) 


1-12 


3-9 


1-10 


4-9 


5-17 


Yield determination (correlated) 


2.5 


2.5 


0.5 


1.0 


3.0 


(uncorrelated) 






1-5 


2-14 


4-9 


PID efficiency 


1-5 


1-5 


6-19 


1-15 


5-9 


Tracking efficiency 


6 


10 


9 


9 


9 


Branching fraction 


1.3 


1.5 


2.1 


5.8 


26.0 


Luminosity 


3.5 


3.5 


3.5 


3.5 


3.5 



4 Cross-section measurements 

The signal yields determined from the data allow us to measure the differential cross- 
sections as functions of pt and y in the range < pt < 8 GeV/c and 2.0 < y < 4.5. The 
differential cross-section for producing hadron species H c or its charge conjugate in bin i, 
d<7j(if c )/dpr ; integrated over the y range of the bin is calculated with the relation 

da,(# c ) 1 N t (H c ^ f + c.c.) 



dp T A PT e i>tot (H c ^ f) ■ B(H C ^ f) ■ C 



int 



where Ap T is the width in p T of bin i, typically lGeV/c, Ni(H c — > f + c.c.) is 
the measured yield of H c and their charge conjugate decays in bin i, B(H C ^- f) 
is the branching fraction of the decay, Ei^ot(H c ^- f) is the total efficiency for ob- 
serving the signal decay in bin i, and C mt = 15.0 ± 0.5nb _1 is the integrated 
luminosity of the sample. The following branching fractions from Ref. [23 are 
used: B(D + ^ K-n + TT + ) = (9.13 ± 0.19)%, B(LT+ -> D°(K-it + )it + ) = (2.63 ± 0.04)%, 
B{A+ -»■ pl<-n + ) = (5.0 ± 1.3)%, and B((D° + D°) K~ti + ) = (3.89 ± 0.05)%, where 
the last is the sum of Cabibbo-favoured and doubly Cabibbo-suppressed branching fractions. 
For the measurement we use the branching fraction of — > K~K + tc + in a ±20 MeV/c 
window around the 0(1020) mass: B(D+ ->■ <j)(K- K + )ti + ) = (2.24 ± 0.13)% [25]. The mea- 
sured differential cross-sections are tabulated in the appendix. Bins with a sample size 
insufficient to produce a measurement with a total relative uncertainty of less than 50% 
are discarded. 

Theoretical expectations for the production cross-sections of charmed hadrons have 
been calculated by Kniehl et al. using the GMVFNS scheme [l}j6] and Cacciari et at, 
using the FONLL approach (7 10 . Both groups have provided differential cross-sections 



as functions of and integrated over bins in y. 
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The FONLL calculations use the CTEQ 6.6 26 parameterisation of the parton den- 
sities. They include estimates of theoretical uncertainties due to the charm quark mass 
and the renormalisation and factorisation scales. However, we display only the central 
values in Figs. |4]-[5j The theoretical calculations assume unit transition probabilities from 
a primary charm quark to the exclusive hadron state. The actual transition probabilities 
that we use to convert the predictions to measurable cross-sections are those quoted by 
Ref. [27] , based on measurements from e + e~ colliders close to the T(AS) resonance: 
f(c-+D°) = 0.565 ±0.032, f(c D+) = 0.246 ± 0.020, f(c -»■ D*+) = 0.224 ± 0.028, 
f (c D+) = 0.080 ± 0.017, and f(c -»■ A+) = 0.094 ± 0.035. Note that the transition 
probabilities do not sum up to unity, since, e.g., f(c — > D°) has an overlapping contri- 
bution from /(c — > D* + ). No dedicated calculation for production is available. The 
respective prediction was obtained by scaling the kinematically similar D* + prediction by 
the ratio f(c D+)/f(c D*+). 

The GMVFNS calculations include theoretical predictions for all hadrons studied in our 
analysis. Results were provided for p^ > 3GeV/c. The uncertainties from scale variations 
were determined only for the case of D° production. The relative sizes of the uncertainties 
for the other hadron species are assumed to be the same as those for the D°. Here the 
CTEQ 6.5 28 set of parton densities was used. Predictions for D° mesons were also 
provided using the CTEQ 6.5c2 29 parton densities with intrinsic charm. As shown in 
Fig. 4a in the phase space region of the present measurement the effect of intrinsic charm 
is predicted to be small. The GMVFNS theoretical framework includes the convolution 
with fragmentation functions describing the transition c — > H c that are normalised to the 
respective total transition probabilities [2J. The fragmentation functions are results of a fit 
to production measurements at e + e~ colliders, where no attempt was made in the fit to 
separate direct production and feed-down from higher resonances. 

To compare the theoretical calculations to our measurements, the theoretical differential 
cross-sections were integrated over the bins and then divided by the bin width Ap^. 
The integration was performed numerically with a third-order spline interpolation of the 
differential cross-sections. 

The measured cross-sections compared to the theoretical predictions are shown in 
Figs. |4]-[5j For better visibility, theoretical predictions are displayed as smooth curves such 
that the value at the bin centre corresponds to the differential cross-section calculated 
in that bin. The data points with their uncertainties, which are always drawn at the 
bin centre, thus can be directly compared with theory. The predictions agree well with 
our measurements, generally bracketing the observed values between the FONLL and 
GMVFNS calculations. 



5 Production ratios and integrated cross-sections 

Charmed hadron production ratios and total cross-sections are determined for the kinematic 
range < p^ < 8 GeV/c and 2.0 < y < 4.5. Bins where the relative uncertainty on the 
yield exceeds 50% (left blank in Tables I5HTU1 of the appendix) are not used. Instead, 
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Figure 4: Differential cross-sections for (a) D°, (b) D + , (c 



and 



(d) Df meson production 



compared to theoretical predictions. The cross-sections for different y regions are shown as 
functions of pp- The y ranges are shown as separate curves and associated sets of points scaled 
by factors 10 -m , where the exponent m is shown on the plot with the y range. The error bars 
associated with the data points show the sum in quadrature of the statistical and total systematic 
uncertainty. The shaded regions show the range of theoretical uncertainties for the GMVFNS 
prediction. 
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Figure 5: Differential cross-sections for A£ baryon production compared to the theoretical 
prediction from the GMVFNS scheme. The error bars associated with the data points show the 
sum in quadrature of the statistical and total systematic uncertainty. The shaded region shows 
the range of theoretical uncertainty for the theoretical prediction. 
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Table 2: Open charm production cross-sections in the kinematic range < < 8GeV/c and 
2.0 < y < 4.5. The computation of the extrapolation factors is described in the text. The first 
uncertainty is statistical, the second is systematic, and the third is the contribution from the 
extrapolation factor. 



Extrapolation factor Cross-section (fib) 



D° 


1.003 


±0.001 


1661 ± 16 ± 128 ± 2 


D + 


1.067 


±0.013 


645 ± 1 1 ± 72 ± 8 


D *+ 


1.340 


±0.037 


677±26± 77 ±19 


Df 


1.330 


±0.056 


197±14± 26 ± 8 


K 


1.311 


±0.077 


233±26± 71 ±14 



Table 3: Correlation matrix of the uncertainties of the integrated open charm production 
cross-sections in the kinematic range < pr < 8GeV/c and 2.0 < y < 4.5. The first column 
restates measured values of the integrated cross-sections. 











a(D*+) 






= 1661 ± 129 fib 










a(D+) 


= 645 ± 74 fib 


0.76 








a{D*+) 


= 677 ± 83 fib 


0.77 


0.73 






°(Pt) 


= 197 ± 31 nb 


0.55 


0.52 


0.53 




a{At) 


= 233 ± 77 fib 


0.26 


0.25 


0.25 


0.18 



the cross-sections are extrapolated from the remaining bins with predictions obtained 
from Pythia 6.4. The extrapolation factors are computed as the ratios of the predicted 
cross-sections integrated over < < 8 GeV/c and 2.0 < y < 4.5 to the predicted cross- 
sections integrated over the well measured bins for each of four tunes of Pythia 6.4: 



LHCb-tune 17 , Perugia 0, Perugia NOCR, and Perugia 2010 [30]. The mean of these 
four ratios is used as a multiplicative factor to extrapolate the sum of the well measured 
bins to the full kinematic range under study. The root mean square of the four ratios is 
taken as a systematic uncertainty associated with the extrapolation. We confirm that this 
procedure gives uncertainties of appropriate size by examining the variance of the ratios 
for individual well measured bins. The resulting integrated cross-sections for each hadron 
species are given in Table [2j 

Accounting for the correlations among the sources of systematic uncertainty, we 
obtain the correlation matrix for the total uncertainties of the integrated cross-section 
measurements shown in Table [3] The ratios of the production cross-sections in the 
kinematic range < px < 8GeV/c and 2.0 < y < 4.5 are given in Table |1J 

Finally, we determine the total charm cross-section contributing to charmed hadron 
production inside the acceptance of this study, < pr < 8 GeV/c and 2.0 < y < 4.5. Com- 
bining our measurements <j(H c ) with the corresponding fragmentation functions f(c — > H c ) 
from Ref. p7] gives five estimates of a(cc) = a(H c )/(2f(c — > H c )). The factor of 2 appears 
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Table 4: Cross-section ratios for open charm production in the kinematic range < px < 8 GeV/c 
and 2.0 < y < 4.5. The numbers in the table are the ratios of the respective row/column. 





a(D°) 


a(D+) 


a{D*+) 




*(£>+) 


0.389 ±0.029 








a(D* + ) 


0.407 ±0.033 


1.049 ±0.092 








0.119 ±0.016 


0.305 ± 0.042 


0.291 ±0.041 




a(At) 


0.140 ±0.045 


0.361 ±0.116 


0.344 ±0.111 


1.183 ±0.402 



in the denominator because we have defined a(H c ) to be the cross-section to produce either 
H c or its charge conjugate. A combination of all five measurements taking correlations 
into account gives 

a(cc) PT<8 Gcv/c, 2.o<j/<4.5 = 1419 ± 12 (stat) ± 116 (syst) ± 65 (frag) \ib, 
The final uncertainty is that due to the fragmentation functions. 

6 Summary 

A measurement of charm production in pp collisions at a centre-of-mass energy of 7 TeV 
has been performed with the LHCb detector, based on an integrated luminosity of 
£i n t = 15 nb" 1 . Cross-section measurements with total uncertainties below 20% have been 
achieved. The shape and absolute normalisation of the differential cross-sections for D°/D°, 
D ± , D*^, Df, and Af hadrons are found to be in agreement with theoretical predictions. 
The ratios of the production cross-sections for the five species under study have been mea- 
sured. The cc cross-section for producing a charmed hadron in the range < p? < 8 GeV/c 
and 2.0 < y < 4.5 is found to be 1419 ± 12 (stat) ± 116 (syst) ± 65 (frag) yib. 
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Appendix 



Measured open charm cross-sections 

Table [5] shows the production cross-sections for A£ baryons integrated over 
2 < pi < 8GeV/c and over the rapidity range of the y bins. The differential produc- 
tion cross-section values (integrated over the y range of the respective bin) plotted in 
Figs. |4]-[5] are given in Tables IBTTTOl 



Table 5: Bin-integrated production cross-sections in (xb for prompt A~£ + c.c. baryons in bins of 
y integrated over the range 2 < p^ < 8 GeV/ c. The first uncertainty is statistical, and the second 
is the total systematic. 



Pt 




y 




(GeV/c) 


(2.0,2.5) 


(2.5,3.0) (3.0,3.5) 


(3.5,4.0) 


(2,8) 


21.4 ±8.1 ±7.2 


49.9 ± 11.6 ± 15.6 62.9 ±7.0 ±18.8 


44.2 ±8.6 ± 13.2 



Table 6: Differential production cross-sections, dr/dpr, in (xb/(GeV/c) for prompt A£ + c.c. 
baryons in bins of pt integrated over the rapidity range 2.0 < y < 4.5. The first uncertainty is 
statistical, and the second is the total systematic. 



Pt 


y 


(GeV/c) 


(2.0,4.5) 


(2,3) 


89.6 ±17.8 ±32.6 


(3,4) 


49.8 ± 7.9 ±15.3 


(4,5) 


22.5 ± 3.1 ± 6.9 


(5,6) 


8.5 ± 1.4± 2.6 


(6,7) 


4.9 ± 0.9 ± 1.5 


(7,8) 


2.4 ± 0.6 ± 0.8 
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